Topological Entanglement Entropy of Z 2 Spin liquids and Lattice Laughlin states 



O 

(N 

Hi 



-f— > 

CO 

03 



i 

c 

o 
o 



> 
in 

o 
o 

o 



Yi Zhang, Tarun Grover and Ashvin Vishwanath 1 
1 Department of Physics, University of California, Berkeley, CA 94720 

We study entanglement properties of candidate wave-functions for SU (2) symmetric gapped spin 
liquids and Laughlin states. These wave-functions are obtained by the Gutzwiller projection tech- 
nique. Using Topological Entanglement Entropy 7 as a tool, we establish topological order in chiral 
spin liquid and Z2 spin liquid wave-functions, as well as a lattice version of the Laughlin state. 
Our results agree very well with the field theoretic result 7 = log D where D is the total quantum 
dimension of the phase. All calculations are done using a Monte Carlo technique on a 12 x 12 lattice 
enabling us to extract 7 with small finite size effects. For a chiral spin liquid wave-function, the 
calculated value is within 4% of the ideal value. We also find good agreement for a lattice version 
of the Laughlin v = 1/3 phase with the expected 7 = log \/3- 

PACS numbers: 



I. INTRODUCTION 

Quantum Spin Liquids(SLs) are states that arise from 
the collective behavior of spins, but are not character- 
ized by a Landau order parameter. They are associated 
with remarkable phenomena such as fractional quantum 
numbers*, transmutation of statistics (eg. fcrmions ap- 
pearing in a purely bosonic model)*^, and enabling oth- 
erwise impossible quantum phase transitions*, to name 
a few. SLs may be gapless or gapped. While current 
experimental candidates for SLs appear to have gapless 
excitations*, gapped SLs arc indicated in numerical stud- 
ies on the Kagome* and honeycomb lattice 7 -. Gapped SLs 
are characterized by topological order - i.e. ground state 
degeneracy that depends on the topology of the underly- 
ing spaced. 

Recently, a novel characterization of gapped SLs has 
emerged using quantum entanglement in terms of the 
Topological Entanglement Entropy (TEE)2r— . This 
quantity takes a fixed value 7 in a topologically ordered 
phase and remarkably can be calculated just knowing the 
ground state wave-function. The entanglement entropy 
of a two dimensional disc shaped region A in a gapped 
phase obeys S2 = clIa — 7, where a smooth boundary 
of length I a is assumed to surround the region. By care- 
fully subtracting off the leading dependence, the constant 
7 can be isolated. It is argued to be a characteristic of 
the phase, 7 = log£>, where D is the quantum dimen- 
sion of the phas o 10 ' 11 . For the abelian states discussed 
here, D 2 is identical to the ground state degeneracy on 
the torus. 

Gapped SLs can be viewed as a state where each spin 
forms a singlet with a near neighbor, but the arrangement 
of singlets fluctuates quantum mechanically so it is a liq- 
uid of singlets. Theoretical models of this singlet liquid 
fall roughly into two categories. In the first, the singlets 
are represented as microscopic variables as in quantum 
dimer and related models* 2 - - — , and are suggested by large 
N calculations* 7 -. Topological order can then be estab- 
lished by a variety of techniques including exact solution 
and most recently quantum entanglement^* - — . In con- 



trast there has been less progress establishing topological 
order in the second category, which are SU(2) symmet- 
ric spin systems where valence bonds are emergent de- 
grees of freedom. Anderson* proposed constructing an 
SU(2) symmetric SL wave-function by starting with a 
BCS state, derived from the mean field Hamiltonian: 
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and Gutzwiller projecting it so that there is exactly one 
fcrmion per site, hence a spin wave-function. Variants of 
these are known to be good variational ground states for 
local spin Hamiltonians (see e.gi**) and are more viable 
descriptions of most experimental and SU(2) symmetric 
liquids. Approximate analytical treatments of projection, 
that include small fluctuations about the above mean 
field state, indicate that at least two kinds of gapped 
SLs can arise: chiral SLs** and Z 2 SL a 17 ' 25 ' 26 . However, 
given the drastic nature of projection, it is unclear if the 
actual wave-functions obtained from this procedure are 
in the same phases. In this paper we use TEE to estab- 
lish topological order for of SU (2) symmetric chiral and 
Z 2 SL wave- functions. We show that the recently devel- 
oped Monte Carlo technique used to study entanglement 
properties of gapless SLs* 7 - can be applied here as well to 
extract TEE for system sizes large enough (144 spins) so 
that it approaches its quantized value. Instead of using 
the more standard Von Neumann entropy, we focus on 
the Renyi entropy, which carries the same contribution 
from the TEE for both non-chiral*£ and chiral^* states 
for a topologically trivial bipartition. The fact that the 
wave-function is a determinant or product of determi- 
nants in these cases allows for its efficient evaluation. 
For a model of of the chiral SL, the calculated TEE are 
remarkably accurate, within few percent of the expected 
log a/2 value. To our knowledge, this is the first clear 
demonstration of topological order via TEE, in SU(2) 
symmetric spin wave-functions. 

We also study lattice versions of the Laughlin v = 1/3 
state, which are obtained by a similar projective con- 
struction, although these are fermionic, not spin wave- 
functions. Again we can extract TEE which is within 
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7% of the expected value to confirm these are in the 
same phase as the Laughlin state, although they differ 
significantly in microscopic structure. 

We note earlier numerical work extracting TEE include 
exact digitalization studies on small systems, looking at 
quantum Hall Laughlin states^ and perturbed Kitaev 
toric code models^. Recently, a quantum Monte Carlo 
study-2£ used TEE to detect Z 2 topological order. In 
contrast to the states studied here, this was a positive 
definite wave-function, with U(l) rather than SU(2) spin 
symmetry. Our wave-function-only approach is ready- 
made for searching for topological order when one has a 
good variational ansatz for a ground state, irrespective of 
whether the state is positive definite or not. Finally, we 
note that Refill studied topological order in 'nodal' Z 2 
SLs by constructing four orthogonal low energy states on 
the torus, and Ref.— studied TEE for the Kitaev model. 

The format of the paper and main results are summa- 
rized in the section below. 

In SectionlDthc TEE 7 is defined, and an algorithm to 
calculate it numerically utilizing the Renyi entanglement 
entropy S 2 is outlined. This is then applied to a series of 
topological phases, the results of which are summarized 
in table HI 

(i) The first is a Chern insulator, built out of a square 
lattice tight binding model at half filling, in which the 
filled band has unit Chern number. For a lattice with 2N 
sites, this is an N body Slater determinant $(ri, . . . , rpj). 
Since this is an integer quantum Hall state, it is not ex- 
pected to possess topological order. Indeed, calculation 
is consistent with a vanishing TEE, see Table Q] first row. 

(ii) The chiral SL wave-function is obtained from the 
wave- function of 2N spinful electrons with this tight bind- 
ing band structure, by projecting out all double oc- 
cupancies, and studied in Section IIII1 The chiral SL 
wave-function can be written as a product of two Slater 
determinants, i.e.\E'(ri, r 2 . . . , rjv) = M^ 2 (r 1: . . . , rjv), 
where M. an unimportant Marshall sign factor. If we 
view up spin as a hardcore boson, then this is the wave- 
function analogous to half filled Landau level v = 1/2 
of bosons. It is therefore expected to have 7 = log 
Indeed, for a particular choice of parameters with a large 
gap, numerical calculation (second and third rows of ta- 
bic HI with different linear dimensions La of the small- 
est subregions involved) yields a value very close to this. 
Detailed finite size analysis obtained by varying the cor- 
relation length of the chiral SL is presented in Section 
IIII1 providing further evidence for convergence to the ex- 
pected value. We note that this wave-function is SU(2) 
symmetric, and non-positivc-definite, since the ground 
state is not time reversal symmetric. 

(iii) Note, the construction of the chiral SL above 
is similar to the Laughlin construction of fractional 
quantum Hall states by taking products of the inte- 
ger quantum Hall states. Extending the construction 
above, one can write wave-functions for N fermions 
^ , 1 / 3 (r 1 , r 2 ■ ■ ■ , r N ) = < f )3 (r 1 , . . . , r N ), a lattice version 
of the Laughlin v = 1/3 state. The entanglement en- 
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log 
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0.99 ± 0.12 


Lattice v = 1/3 


log 
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1.07± 0.05 


Lattice v = 1/4 


log 


\/4 


1.06 ± 0.11 


Z 2 SL L A =4 


1°!- 


-2 


0.84 ± 0.13 



TABLE I: Comparison between calculated TEE and expected 
value from field theory (second column) for topological phases. 
The * denotes that the calculated value is not divided by the 
expected value since the latter vanishes. 



tropy calculation for ^1/3 agrees well with what is ex- 
pected for the topological order for u = 1/3 Laughlin 
state, indicating it is in the same phase, despite not be- 
ing constructed from lowest Landau level states. Note, 
since they differ significantly in microscopic detail from 
the Laughlin state, wave-function overlap is not an op- 
tion in establishing that they are in the same phase. 
Also, calculating entanglement spectral is currently not 
feasible for these wave-functions, thus TEE appears to 
be the ideal characterization. Similarly, the lattice ana- 
log of Laughlin v = 1/4 state for bosons, obtained via 
^1/4 ('"li r 2 ■ ■ ■ , Hv) = ^ 4 0"ii ■ • ■ , tn), is found to have a 
TEE close to the expected 7 = log \/4, as discussed in 
Section [TV] 

(iv) Finally, we construct a fully gapped Z 2 SL wave- 
function on the square lattice. For the largest system 
sizes we considered, the calculated 7 is 84% of the ex- 
pected log 2 value (last row in table [I]). The difference is 
ascribed to larger finite size effects, as discussed in Sec- 
tion |V] 



II. TOPOLOGICAL ENTANGLEMENT 
ENTROPY AND VARIATIONAL MONTE 
CARLO METHOD 

A. Renyi entropy and topological entanglement 
entropy 

Given a normalized wave- function |$) and a partition 
of the system into subsystems A and B, one can trace out 
the subsystem B to obtain the reduced density matrix on 
A: pa = Ttb 1$) ($|. The Renyi entropies are defined 
as: 

S n = ^— log (Trp n A ) (2) 

1 — 71 

Taking the limit n — > 1, this recovers the definition 
of the usual von Neumann entropy. In this paper we 
will focus on the Renyi entropy with index n = 2: 
S 2 = — log (TV (pa)), which is easier to calculate with 
our Variational Monte Carlo(VMC) method^. 

For a gapped phase in 2D with topological order, a 
contractible region A with smooth boundary of length 
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I A, the Area Law of the Renyi entropy becomes: 
S 2 = al A -~f 

where we have omitted the sub-leading terms. Al- 
though the coefficient a of the leading boundary law' 
term is non-universal, the sub-leading constant 7 is uni- 
versal, and this TEE is a robust property of the phase of 
matter for which |$) is the ground state. It is given 
by 7 = log-D, where D is the total quantum dimen- 
sion of the mode l 10 ' 11 , and offers a partial characteri- 
zation of the underlying topological order. When re- 
gion A has a disc geometry, it has been shown that 7 
for different Renyi indices n are identical for both chiral 
and non-chiral state a 28 ' 29 . A simple limit where this is 
readily observed 9 ^ is in a model wave-function of a Z 2 
SL, which is an equal superposition of loops (Z 2 elec- 
tric field). This is achieved as a ground state in Ki- 
taev's toric code model^i. The Schmidt decomposition 
into wave-functions in regions A and B can be indexed 
by the configuration of electric field lines piercing the 
boundary of the disc. If i = 1, 2, ... I are I bonds going 
through the boundary between region A and B, the pres- 
ence (absence) of electric field lines on bond i is denoted 
by qi = 1 (qi = 0). Since the loops are closed, we require 
Si Qi = even. It can be shown that the wave-function is 
simply an equal weight decomposition indexed by all pos- 
sible configurations of q^. There are C = 2 /_1 of them, 
the global constraint of closed loops accounting for the 
missing factor of 2. Then: 



I*) = 



f 



91 + 



■■+qi even 



91— gi> 



If 



qi—qii 



This implies^ there are C equal eigenvalues of the re- 
gion A density matrix, each equal to 1 jC . The Renyi en- 
tropy from EqnUJis: S„ = log C*^™- 1 ) = (7-1) log 2. 
Thus 7 = log 2 from the definition above, if we identify I 
with the length of the boundary. Note, this follows inde- 
pendent of the Renyi index n, and is the expected value 
for a Z 2 gauge theory with quantum dimension D = 2. 

Practically, it is not convenient to extract the sublead- 
ing constant by fitting the expression above, particularly 
on the lattice where edges frequently occur. Instead, one 
may use a construction due to Levin and WeniS, or Ki- 
taev and PreskilUi, that effectively cancels out the lead- 
ing term and exposes the topological contribution. We 
use the latter, which requires calculating entanglement 
entropy for a triad of non overlapping regions A, B, C, 
and their various unions, and then constructing: 
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Sa + Sb + Sc — Sab — Sac 



S 



BC 



Sabc 



here, any S n can be used, and we choose to use 5*2, 
since it can be easily calculated. This guarantees that the 
contributions of boundaries and corners cancel when the 
dimensions of individual regions A,B,C is much larger 
than the correlation length. 



B. Variational Monte Carlo method for Renyi 
Entropy 

In this section we briefly review the VMC algorithm 
for calculating Renyi entropy S 2 27]35 - Consider the con- 
figurations \a x ) = \a)\b), \a 2 ) = \a')\b'), Ift) = \a')\b), 
IA2) = l a )l^')i where |a) and \a') have their support only 
in the subsystem A while \b) and \b') are in subsystem B. 
Following Refi^, we define an operator SwapA that acts 
on the tensor product of two copies of the system and 
swaps the configurations of the spins belonging to the 
A subsystem in the two copies i.e. SwapA|ai) ® \a 2 ) = 
<8 1/52) ■ The Renyi entropy S 2 for the bipartition A 
and B can be expressed in terms of the expectation value 
of SwapA with respect to the wave- function |$) <£> |<i>): 



Sa = - log (trp 2 A ) = - log (Swap A ) 



(3) 



< SwapA > may be re-expressed as a Monte Carlo 
average: 



(Swap A ) = ^ P ai p a2 f («i,a 2 ) 



(4) 



where the weights p Ql = |(aj|$)| /J2 a , \( a i\®)\ are 
normalized and non-negative while the quantity to be 
averaged over the probability distribution p ai p a2 i s: 



1 (ai '" 2)= < ai |*> <a 2 |*> 



(5) 



Therefore, one can calculate the Renyi entropy using 
VMC method. This technique is particularly suited for 
projected wave-functions since the projection is rather 
easy to implement in a VMC algorithm^. As shown 
in ReL2£ VMC algorithm correctly reproduces the exact 
results for free fermions with an error of less than a few 
percent. 

We further facilitate our calculation with an algorithm 
that we referred to as the sign trick^l. It offers simpli- 
fication and reduces computational cost. Basically, we 
separate (SwapA) as a product of two factors, which may 
be independently calculated within VMC method: 



(Swap A ) = (SwapA.mod) (SwapA.si, 



Q1Q2 



P ai p a2 \f (ai,a 2 ) 



E~ „i<P(a 1 ,a 2 ) 
Hai,a 2 t> 

_OL\OL 2 



The first factor is the Renyi entropy of a sign prob- 
lem free wave-function \4> ai \- The second term is 
the expectation value of the phase factor e l ^ ai ' a2 ' = 
^yy/ii^/ l^a^aa^i^l with probability distri- 
bution p aua2 = 10*^*2 0^0^1/ E l^Ka faith I- 

ot\a 2 

Both factors can be calculated in a more efficient manner 
and most importantly, have much smaller errors than the 
direct calculation of (SwapA). 



4 



III. ENTANGLEMENT ENTROPY FOR A 
CHIRAL SPIN LIQUID 

In this section we calculate the Renyi entropy 5*2 and 
TEE 7 for a chiral SL21 



A. Projected wave-function for chiral spin liquids 

The chiral SL is a spin SU(2) singlet ground state, that 
breaks time reversal and parity symmetr y 24 ! 37 . A wave- 
function in this phase wave function may be obtained us- 
ing the slave-particle formalism by Gutz wilier projecting 
a d + id BCS stated. Alternately, it can be obtained by 
Gutzwillcr projection of a hopping model on the square 
lattice. This model has fermions hopping on the square 
lattice with a it flux through every plaquette and imagi- 
nary hoppings across the square lattice diagonals: 



H 



«ifc» 



t,fi fk 



(6) 



Here i and j are nearest neighbors and the hopping 
amplitude Uj is t along the y direction and alternating 
between t and —t in the x direction from row to row; and 
i and k are second nearest neighbors connected by hop- 
pings along the square lattice diagonals, with amplitude 
Aik = iA along the arrows and A^ = — iA against the 
arrows, see Fig. 1. The unit cell contains two sublat- 
tices A and B. This model leads to a gapped state at 
half filling and the resulting valence band has unit Chern 
number. This hopping model is equivalent to a d + id 
BCS state by an SU (2) Gauge transformation 3 ^. We use 
periodic boundary conditions throughout this section. 

The unprojected ground state wave- function \<p) is ob- 
tained by filling all the states in the valence band (e& < 0) 



|0) where 7* = MWEfLs 



V'-bCOS/tb s e lk ' TB is the creation operator for a valence 

electron with spin s and momentum k, and ^A(k)(tpB(k)) 
is the wave- function on sublattice A(B). The projected 
wave-function that corresponds to the chiral SL is ob- 
tained as |$) = P\<j))i where P is the Gutzwillcr pro- 
jection operator that projects the wave-function to the 
Hilbert space of one electron per site. This is imple- 
mented by restricting \a) to the Hilbert space of spins 
i.e. one particle per site. Due to the fact that this 
Hamiltonian contains only real bipartite hoppings and 
imaginary hoppings between the same sublattices and 
preserves the particle-hole symmetry, this wave-function 
(a|(f>) can be written as a product of two Slater deter- 
minants M.Det(Mij) 1 , where A4 is just an unimportant 
Marshall sign factor, and: 
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FIG. 1: Illustration of a square lattice hopping model con- 
nected with a d+id superconductor. While the nearest neigh- 
bor hopping is along the square edges with amplitude t (— t 
for hopping along dashed lines), the second nearest neigh- 
bor hopping is along the square diagonal (arrows in bold), 
with amplitude +iA (—iA) when hopping direction is along 
(against) the arrow. The two sublattices in the unit cell are 
marked as A and B. 



Tj is the coordinates of the up spins in configuration a, 
and fcj is the momentums in the momentum space. The 
Renyi entropy Si of this wave-function can be calculated 
by VMC method detailed in the last section. 

For an accurate calculation of TEE 7, it is important 
that the sublcading terms in the Eqn. [3] be much smaller 
than the universal constant 7 itself. This finite size error 
is suppressed when the excitation gap is large and cor- 
relation length is shorter than the system typical length 
scale. Note that the mean field gap is given by 8A for 
|A| < 0.5t and 2t v /8 - (t/A) 2 for |A| > 0.5t. To mini- 
mize the finite size effect, we take A = 0.5t unless oth- 
erwise specified, so that the gap is large in both units 
of t and 2A, and our calculation estimates a correlation 
length of £ - 0.45. 



B. 



Establishing Topological Order in chiral SL 
wave- function 



calculate 



the TEE 7 using the 



In this section we 
Kitaev-Preskill scheme 

We study system with total dimensions 12 x 12 lattice 
spacings in both directions with periodic boundary con- 
ditions. We separate the system into La x La squares A 
and B and an La X 2La rectangle C, see Fig [2] For this 
particular geometry, TEE is simply given by: 



Ma- 



{[i>A(h) + Mb)) + (-i) % [Mh) - i>B{h)]}e lk ^ 



7 = 25; 



2. A 



2S- 



2, AC 



~>2,ABC 



(7) 
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FIG. 2: The separation of the system into subsystem A, B, 
C and environment, periodic boundary condition is employed 
in both x and y directions. 
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FIG. 3: Illustration of finite size effect: chiral SL TEE 7 as 
a function of 2A/t proportional to the relative gap size for 
characteristic system length La = 3. The larger the gap, the 
closer the data approaches the ideal value. For comparison, 
TEE 7 for chiral SL at L A = 4 and 2A/t = f.O is shown. On 
the same plot, TEE 7 of a lattice version of v = 1/3 Laughlin 
state at La = 3, 2A/t = 1.0 is also shown. The dashed lines 
are the ideal TEE values of 7 = log (^2) for the chiral SL 
and log (\/3) for the v = 1/3 Laughlin state. 



where we have used the fact that £2, a = 5*2, b, S2,ab = 
S2,c and 5*2, ac = S2,bc owing to the reflection and 
translation symmetry of the wave function. This sim- 
plifies the measurement of TEE into the measurement of 
£2 for only three subsystems A, AC and ABC. 

We use the unprojected wave-function as a benchmark 
for extraction of TEE, which is non-interacting and hence 
exactly solvable. For an La = 3 system, the VMC calcu- 
lation gives 7 = —0.0008 ± 0.0059, in agreement with the 
absence of topological order and correspondingly vanish- 
ing TEE(table[J first row). 

The Gutzwillcr projected wave-function is believed to 
be a chiral SL which can be thought of as a Laughlin 
liquid at filling v = 1/2. Using VMC method, we find 7 = 
0.343±0.012 for an La = 3 system and 7 = 0.344±0.043 
for an La = 4 system, both are in excellent consistency 
with the expectation of 7 = log (V%) = 0.347 for its 
ground states' two fold degeneracy, see table U second 
and third rows and also Fig[3] 

We also want to point out that Gutzwillcr projection 
qualitatively changes the system ground state's topolog- 
ical and quantum behavior from the mean field result. 

On the other hand, by lowering the ratio of 2A/t and 
correspondingly the gap size the correlation length in- 
creases and the finite size effects from sublcading terms 
become more important. See FigJ3] for the approach 
of the extracted TEE 7 to its universal value of 7 = 
log (V%) as we lift the gap size controlled by 2A/t for a 
system with typical length scale La = 3. The finite size 
analysis and the above consistency between La = 3,4 
confirm that finite size effect is small for our chosen sets 
of parameters for the system sizes we study. 



IV. A LATTICE VERSION OF THE LAUGHLIN 
STATE 

Using VMC method, we further study the situations 
where the wave-function is the cube or the fourth power 
of the Slater determinant of the Chern insulator. For 
example, consider the wave-function: 

*1 /3 (JT, r 2 . . . , r N ) = $ 3 (r 1; . . . , r N ) (8) 

where $ is the Chern insulator Slater determinant de- 
fined above. Clearly, the product is a fermionic wave- 
function, since exchanging a pair of particles leads to a 
sign change. This is similar in spirit to constructing the 
corresponding Laughlin liquid of m — 3 of fcrmions, by 
taking the cube of the Slater determinant wave-function 
in the lowest Landau level ip(zi, . . . , zjv) = Ili<j( z i — 

z j) e 41 B )• However, unlike the canonical Laughlin 
state, composed of lowest Landau level states, these are 
rather different lattice wave-functions. An interesting 
question is whether the lowest Landau level structure is 
important in constructing states with the topological or- 
der of the Laughlin state, or whether bands with identical 
Chern number is sufficient, as suggested by field theoretic 
arguments. 

To address this we calculate TEE and compare with 
expectation for the Laughlin phase. Again we choose 
La = 3 in our VMC simulation, and obtain: 7 = 0.5894± 
0.0272 for the m = 3 wave-function, in reasonable agree- 
ment with the ideal value 7 = log (V3) = 0. 549 (table U 
fourth row). 
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We also considered the fourth power of the Chern in- 
sulator Slater determinant: 

* 1/4 (ri, r 2 . . . , r N ) = $ 4 (n, ...,r N ) (9) 

this is a bosonic wave-function, that is expected to 
be in the same phase as v = 1/4 bosons. Indeed we 
find with La — 3 in our VMC simulation, 7 = 0.732 ± 
0.076, consistent with ideal value that must be realized in 
the thermodynamic limit of this phase: 7 = log (\/4) = 
0.693(tableU fifth row). 

These results offered direct support for the TEE for- 
mula 7 = log I? as well as their validity as topological 
ground state wave-functions carrying fractional charge 
and statistics. The lattice fractional Quantum Hall wave- 
functions discussed here may be relevant to the recent 
studies of flat band Hamiltonians with fractional quan- 
tum Hall states^l 4 ^ 



V. ENTANGLEMENT ENTROPY OF A Z 2 SPIN 
LIQUID 

With the projected wave-function ansatz, we may also 
construct a topological Z2 SL by projecting another 
mean-field BCS state, given by the specific BdG Hamil- 
tonian on a square lattice as the following^: 

H = -JL (ft ^3 + h - c ) + Yfi ^^* 

where ipi = f/t)/|) • t 1 ' 2 ' 3 are Pauli matrices. The 
second term is related to chemical potentials, we set 
a 2 ,' 3 = 0, with aj fixed by the conditions /^V 1 ' 2 ' 3 ^) = 0. 
Matrices ^ connecting nearest and next nearest neigh- 
bors: 

3 

^i^i+x l^i^i+y T 

Hi,i+x+y = + 
Hi,i-x+y = TIT 1 — At 2 

This mean field model is readily solvable, with disper- 
sion: 

E k = y^ + lAH 

e k = 2(cos(fc x ) + cos(k y )) 
Afc = 2-q [cos (k x + k y ) + cos (k x — k y )\ + Oq 
— 2iX [cos {k x + k y ) — cos (k x ~ k y )] 

We choose 77 = A = 1.5 for a large gap and our calcu- 
lation estimates that the correlation length is as short as 
£ ~ 1.3 lattice spacings. 

The VMC algorithm need little change^, except 
that instead of Slater determinants product, the wave- 
function for spin product configuration \a) is given by 



<j) a = (a|$) = det(ay) 

here a%j = a (r^ — Tj i) is the Fourier transform of the 
superconducting pairing function/^, r^-j. and Tj^ are the 
coordinates of the up-spins and down-spins, respectively: 

Ik \E k + e k \ 

For numerical simulations we again study 12 x 12 lattice 
spacing systems and separate the system into subsystems 
including La x La squares A and B and La x 2La rect- 
angle C, again see Fig[2j The TEE 7 is given by EqJ7| 
as before. First, we use the unprojected BCS state as 
a benchmark, for which we expect a result of 7 = 0.003 
from an exact solution (since the unprojected state is 
a free particle ground state, one may use the correla- 
tion matrix method^) and consistent with its absence of 
topological order. Indeed, we obtain 7 = 0.012 ± 0.062 
using VMC method, the almost vanishing value of 7 is 
consistent with the expected value, which also serves as 
a check on our Monte Carlo calculations. 

On the other hand, the projection qualitatively al- 
ters the topological properties of the system, and for 
simulation accuracy and efficiency, we employ the 'sign 
trick' from Ref.— . For an La = 4 system, the VMC 
calculation gives 7 = 0.584 ± 0.089. This is roughly 
consistent the Z2 SL which has D 2 = 4 sectors and 
7 = log (D) = log (2) ~ 0.693. The TEE is found to 
be about 84% of the expected value(table H] last row). 
Other studies on Z2 phases eg. the quantum Monte Carlo 
on a Bosc Hubbard model in Refi^ have also found val- 
ues that underestimate the topological entropy (75% of 
the expected zero temperature value in that case). The 
smaller than expected value is probably due to spinon ex- 
citations with a finite gap, causing breaking of Zi electric 
field lines over the finite system size La = 4 we consider. 
Indeed, spin correlations decay more slowly for the Zi 
state, as compared to the chiral SL, which also arrives 
closer to its expected 7 value. Consistent with this fact 
is the observation that for a smaller system size La = 3 
where the finite system size has a larger impact, VMC 
calculation leads to a value of 7 = 0.446 ±0.119, which 
is further away from the ideal value. 

VI. CONCLUSION 

In this paper we studied entanglement properties of 
candidate wave-functions for SU(2) symmetric gapped 
SLs and Laughlin states, and established their topo- 
logical order using the notion of TEE. We studied two 
classes of SLs: 1) Wave-functions that describe quantum 
Hall states and are obtained from the wave-function of a 
Chern insulator by taking multiple copies of it 2) A Z2 
SL state that is obtained by Gutzwillcr projecting a fully- 
gapped BCS superconductor. These wave-functions have 
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long been used as ansatz for exploring SLs states and it is 
reassuring that topologically ordered states can be good 
variational ground states for realistic Hamiltonians. Our 
method is directly applicable to any wave-function that 
can be dealt within VMC method and would be espe- 
cially useful in cases where one is dealing with a Hamil- 
tonian that has Monte Carlo sign-problem and only has 
a variational ansatz for the corresponding ground state. 
We also note that since the quantum Hall wave-functions 
we study are not constructed from the lowest Landau 
level but rather from the band structures that have non- 
zero Chern number, our results are also relevant to the 
recently discovered quantum Hall physics in flat band 
Hamiltonian a 39 ' 40 . 

Let us consider a few problems where our method may 
find immediate application. Firstly, it would be inter- 
esting to apply our method to Z-i SLs that have gaplcss 
nodal spinons. These SLs are obtained by Gutzwiller 
projecting a nodal BCS state. We note that in this case 
one would find an additional contribution to the sublead- 
ing constant part of the entanglement entropy that comes 
from the gaplcss spinons. Though we believe it might still 
be possible to separate the total contribution of the con- 



stant term into a topological constant and a term that 
comes from the gapless spinons only. It would also be 
interesting to study wave-functions that are expected to 
have non-abelian quasi-particles such as SU(2)k quantum 
Hall wave-functions^. Thirdly, since VMC techniques 
can be used for wave-functions defined in the continuum 
as well, it might be interesting to study TEE of quantum 
Hall wave-functions (and their descendants such as time- 
reversal invariant fractionalized topological insulators^, 
) defined directly in the continuum. 

Finally, we note that one limitation of our method is 
that it can't be used to calculate TEE for SLs where 
the gauge fields arc coupled to bosonic (rather than 
fcrmionic) spinons. This is because VMC techniques are 
not very efficient when dealing with wave-functions that 
are written as pcrmanents (in contrast to determinants). 
It would be interesting to see if the recent VMC calcula- 
tion for a SU(2) symmetric bosonic Sl4^ can be pushed 
to bigger system sizes so as to establish topological order 
in such wave-functions. 
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